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We theoretically show that bright six-partite continuous- variable entanglement can be generated 
D ' using cascaded four-wave mixing effects of third-order nonlinearity atomic systems above thresh- 

old. The six-partite continuous-variable entanglement among the six cavity fields with different 
frequencies is analyzed by applying optimized inseparability criteria proposed by Van Loock and 
Furusawa. 
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I. INTRODUCTION 

Quantum entanglement is one of the most mysterious quantum phenomena and has become an indispensable 

■ resource for quantum communication and quantum computation. In the past few years, there has been much effort 
' devoted to the study of continuous-variable (CV) quantum communication , of which the basis is the multi- 

computation and controlled quantum teleportation. It was predicted theretically and demonstrated experimentally 

i— I 

that multi-color CV entanglement can be produced by nondegencrate optical parametric oscillator (OPO), optical 
supcrlattice and second-order nonlinearity media [3|-|6|. Alternatively, CV entanglement can also be produced in 

■ coherent atomic systems, whose entangled beams have narrower linewidths and longer coherence times; in particular, 
• they can be used for free-space quantum teleportation 0-S], such as, from a satellite to a ground station. Such 



entanglement using four- wave mixing (FWM) has been predicted by Ref. [10h12| |. All those results above endow 
us the feasibility and necessity to realize our scheme. Based on those processes, we extend the scheme to generate 
six-partite CV entanglement by intracavity four-wave-mixing cascaded with double such generations in a four-level 
atomic system. And we theoretically demonstrate such possibility of multy-partitc entanglement using third-order 
nonlinearity according to the inseparability criterion for multipartite CV entanglement proposed by van Loock and 
Furusawa. 

This paper is arranged as follows. The physical model and its hamiltonian are discussed in Sec. II. Sec. Ill gives the 
equations of motion of the system and Sec. IV provides a linearized fluctuation analysis to calculate the measurable 
output fluctuation spectra. These output spectra are further discussed in Sec. V and they demonstrate the violation 
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FIG. 1: (Color online) (a) Experimental setup for the proposed scheme, (b) Relevant energy level diagram of the four-level 
atomic system. This involves three FWM processes, (c) Frequencies of the pumps and cavity fields. All the processes and 
photons are marked with different colors for simplicity. 



of the optimized VLF criteria and some entanglement characteristics. At last, we draw some brief conclusions in Sec. 
VI. 



II. PHYSICAL MODEL AND ITS HAMILTONIAN 

We consider an ensemble of independent four-level atoms inside a resonant optical cavity. The proposed experimental 
setup for this scheme is shown in Fig. 1(a). A vapor cell containing a four- level atomic system is placed inside a 
Fabry P'erot cavity, using lasers whose frequency are u) p i and uj P 2 as pumps, where the virtual energy level diagram 
of three cascaded FWM effects are shown in Fig. 1(b). Those generated six beams are transmitted out through the 
coupling mirror 2 with partial tramittion coefficient, then being splitted by the prism. Spatially separated beams can 
be measured using method as Ref. [tH to determine the entanglement relationship. 

As shown in Fig. 1, the lower levels, |1), |2) and |3), of the atoms are coupled to the excited level |4) by two 
driving lasers of frequencies lu p i, u P 2. There are three FWM processes which are indicated in Fig. 1(b), that 1 f and 
1 l compose the first FWM process, where 1 f implies io p \ and ui P 2 stimulated from |1) and |2) to virtual levels, 1 4- 
implies u> s i and falling from virtual levels to |2) and |1), respectively. Similar are to 2 f, 2 J,, 3 t an d 3 \.. And we 
suppose all the virtual levels are detuned with |4). 

In present scheme, we only concentrate on the entanglement characteristics among optical fields instead of these 
between atoms and optical fields. So the interaction Hamiltonian of the three cascaded FWM processes could be 
written as 

H m t = i?i[fcia p iap2aiiali + k 2 a sl a pl al 2 al 2 + haaa^a^a^} + H.c. , (1) 
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where ki (i=l,2,3) are the dimensionless nonlinear coupling coefficients and are commonly taken to be a real number, 
aj and a\ are corresponding to the light field generation and annihilation operators repectively. 
Considering the driving fields in the cavity, the Hamiltonian for the pump beams is given by 

Hp Ump = (hleta'p! + e 2 a£ 2 ] + H.c. , (2) 

where ej (i=l,2) is the classical pump and coupling laser amplitude. The loss of the zth mode [ijjl in the cavity can 
be written as 

Lip = ji[2aipal - ajaip - pat a*] , (3) 

where ji stand for the damping rates for the corresponding cavity modes which are related to the amplitude trans- 
mission coefficients. 



III. EQUATIONS OF MOTION FOR THE FULL HAMILTONIAN 



We now consider the full physical system, where the four-level atomic system is contained inside a pumped resonant 
Fabry-P'erot cavity. The master equation for the six cavity modes is 

6 
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The equations of motion for the six cavity modes can be obtained by solving the Fokkcr-Planck equation in the 
P-representation [3] as 
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where a can be expressed as [a P 2, a p i, an, a s x, a^, «s2, o;* 2 , a pl , a* 1; a* 2 , a* 2 ] and rj be expressed as 

[771 (i), 772 (t)...77n(t), r]i2(t)] T , where elements in the vector 77 are real noise terms and have (rji(t)) = 0, (rji(t)rjj (£')) = 
SijSit-t'). 

And F is the main part of system's evolution, as the form of F = [f , f*], where 



\ 



(6) 



£2 - l P 2Ct P 2 - kia. vl a s \Q4\ - k 2 a* 1 a pl a l2 + k 3 a* s2 a sl a pl 
ei - l P ia P i - kia* 2 a s iaii - k 3 a* sl a P 2a S 2 + k 2 a* 2 a l ia P 2 
-Jaan + kia* sl a p ia p2 - k 2 a* 2 a p iai2 
— "f s ia s i + kia* 1 a p ia P 2 — k 3 a pl a P 2a s2 
-H20ti2 + k 2 a* pl a P 2an 
\ -7s2C*s2 + k 3 a* 2 a p ia s i J 

B contains the coefficients of the noise terms and is not calculated here. B can be obtained by factoring the diffusion 
matrix D according to D = BB T . 
We express D as 



D 




(7) 
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where is a 6x6 zero matrix and the nonzero block is given by 
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IV. STABILITY ANALYSIS 



Now we conduct the stability analysis of the system in order to analyze the output spectral correlations. By 
neglecting the noise terms in Eq. (5), we obtain a set of classical equations for the mean values: 
and from these we can obtain steady-state solutions. 

a 



at 
d_ 

~dt 



a p2 — £2 - 7p2" P 2 - kia^agian - k 2 a* 1 a p iai2 + k 3 a* 2 a s ia p i , 



Opi = ex — Jpioipi - kia* 2 a s ian - k 3 a* 1 a p2 a s2 + k 2 a* 2 a>na p2 , 



d_ 

dt 
d 

dl ail 



-liiun+kia* sl a p ia p2 - k 2 a* 2 a pl a i2 , 
-7 s ia s i+fc 1 a* 1 a pl a P 2 - k 3 a* pl a p2 a s2 , 



(9) 



dl as2 
d 



-lao-ii + k 2 a pl a P 2aii 



gj. a i2 = -1siol s2 + k 3 a* 2 a pl a sl , 

In order to simplify the calculation, the same damping rate and the same classical pumping laser amplitudes in 
the cavity are assumed (i.e., 7 P i = 7 p2 = j a , jn = j s i = 76, 7 S 2 = 7i2 = 7c and e% = e 2 = e ). Similarly, the same 
nonlinear coupling coefficient are assumed (i.e., k 2 = k 3 ). 



Threshold can be obtained if k\ > 2k 2 yj^b/lc- The values of the two threshold are expressed as following respecively, 

(10) 



e th =la\l —2 and e th = 7a I/ 9T2 
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When eth < e < e th , there is only one set of nonzero solutions for all six modes. If e th < e, there will be two sets of 
nonzero solutions for all six modes. When there exists a threshold the analytical expressions of stationary solutions 
for different values of the pump amplitude are 

(i) e < eth 
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(iii) e > e 
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A = A a = ^ or A t =A a = ^,(i = pl, P 2); 

la la (13) 



If fci < 2k2y / ib/lc, there is no threshold in this system. On this condition the signal and the idler modes will 
not be excited no matter how large the pump amplitude is. Thus, when a threshold docs not exist the analytical 
expressions of the stationary solutions are 

At = — , (i = pl,p2); Ai = 0,(i = il, si, »2, s2) . (14) 

la 

In the following, we decompose the system variables on into their steady-state values A; and small fluctuations 
around the steady-state values bon as on = Ai + dai. Then equation (5) can be linearized as 

JUa = -MSa + Br] , (15) 

where 5a = [<5a p2 , 6a p i, San, Sa s i, 6ati2, 8a S 2, 8a* 2 , 5a* ± , 5a* 11 8a* sll 6a* 2 , Sa* 2 ] T ; M, the drift matrix, is given by 
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If the requirement that the real part of the eigenvalues stay positive is satisfied, the fluctuation equations will 
describe an Ornstein-Uhlenbeck process for which the intracavity spectral correlation matrix is 



-k 2 A a A b 







S(w) = (M + iwI)- 1 BB T (M T - iuiy 1 



(19) 
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All the correlations required to study the measurable extracavity spectra are contained in this intracavity spectral 
matrix. The output fields can be obtained along with the standard input-output relations 
output fields spectral variances and covariances have the general form 



17| . In particular, the 



STUM =2Sms Xt , Xj (u) 



(20) 



Similar expressions can be derived for the Y quadratures. 



V. MULTIPARTITE ENTANGLEMENT AND OUTPUT SPECTRA 

In order to investigate multipartite entanglement, and particularly show that the system under consideration demon- 
strates true six-partite entanglement, we define quadrature operators for each mode as 

X, = a t + a\ , Y = -i{cn - at) , (21) 

such that [Xi, Yj\ = 2i. The conditions proposed by van Loock and Furusawa [l^, which are a generalization of the 
conditions for bipartite entanglement, are sufficient to demonstrate multipartite entanglement. 
Using the quadrature definitions, the six-partite inequalities are 

V(X i2 - X p i) + V(g p2 Y p2 + Y p i + gaYa + g sl Y sl + Y i2 + g s2 Y s2 ) > 4 , 
V(X pl + X sl ) + V(g p2 Y p2 + Y pl + gaYa - Y sl + g i2 Y l2 + g s2 Y s2 ) > 4 , 

V{X sl - Xa) + V(g p2 Y p2 + g pl Y pl + Yi + Y sl + g l2 Y l2 + g s2 Y s2 ) > 4 , (22) 
V{X n + X p2 ) + V(-Y p2 + g pl Y pl + Y n + g s iY sl + g l2 Y l2 + g s2 Y s2 ) > 4 , 
V(X p2 - X s2 ) + V(Y p2 + g pl Y pl + gaYa + g sl Y sl + g l2 Y l2 + Y s2 ) > 4 , 

where V(A) = (A 2 ) — (A) 2 denotes the variance and gi are arbitrary real parameters used to optimize the violation 
of these inequalities. It is important to note that in the uncorrelated limit these optimized VLF criteria approach 4. 
Hence, without optimization, some entanglement which is presented may be missed. 

From now on, we numerically calculate the values of VLF inequalities, and these are the quantities that can be 
measured in experiments. Note that in our symmetrical situation, inequalities i2-pl and s2-p2 are equal; the same 
are to pl-sl and il-p2. To simplify, we write the subscript sl-il as A, pl-sl and il-p2 as B, i2-pl and s2-p2 as C 
separately. 

When below the threshold or without threshold, the powers of the pump are much higher than others and the 
steady-state values of the other output beams are all equal to zero . Therefore, the quantum characteristics of the 
pump and other beams in such situation are not considered. 

In Fig 2, we plot the minimum of the inequalities as a function of frequency normalized to j a with ^ a ~ lb ~ 1c ~ 
0.03, k\ = 2k 2 = 1, and e = 1.5e t h- In this case, k\ = 2k 2 \J 7e, /7 c so that the two supposed thresholds coincide to 
be one threshold. It is obvious that the minimal values of inequalities are all less than 4 in a wide range of analysis 
frequencies, which is sufficient to demonstrate that all the six beams are CV entangled with each other. In general, 
the smaller the value of the inequality, the stronger the degree of entanglement and quantum correlation is. For large 
frequencies the inequalities approaches 4, which is the result of our optimization. In Fig. 2, we can also see that the 
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FIG. 2: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7„ with 
fa = 7b = 7c = 0.03, ki = 2k 2 = 1, and e = 1.5e t h- There 
is only one threshold in this case. 



FIG. 3: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7 a with 
7a = 7t = 7c = 0.03, fci = 1, fe = 0.4, and e = 1.2e th . 
Here e th = 2e t h and eth < e < e th . 



value of Va is comparatively smaller than Vb , which indicates that the degree of entanglement of second-class output 
is bigger. It is similar to those situations in the four wave mixing process. 

In Fig. 3, we plot the minimum of the inequalities as a function of frequency normalized to "f a with k% = 0.4, e = 
X.ltth and the other parameters remain unchanged. Again, the minima of the inequalities arc less than 4 when above 
threshold, and it verifies that the all field modes are entangled. 

However, we are unable to decide the real entanglement when e > e th because of the two possible stationary solutions. 
Because the solution of ordinary differential equations (ODE) usually depends on the initial value of variables. While 
the variables of ODE here are set initially as Gaussian states and randomly generated at each specific situation. It 
should be expected that all those variables will fall into one of the stationary solutions as the time approaches infinity, 
with a certain distribution of these two stationary solutions. Therefore, we calculate the minimum of the inequalities 
beyond the larger threshold e th using each stationary solution separately. 

In the following, we plot the minimum of the inequalities as a function of frequency normalized to j a using each 
stationary solution at different e, with other parameters unchanged. We find an interesting result. As shown in Fig. 
4 and Fig. 5, when e is a little bigger than e th , the stationary solution A a = tth/la shows a relatively competitive 
result, just as results in Fig. 2 and Fig. 3. The stationary solution A a = e th /"/ a does not work effective and the 
minimum value of Vc is almost near 4. It is because the output beams other than the pump is so lightly occupied if 
A a = e t ij/7o- In other words, A4 (i=sl,il,s2,i2) is close to zero. 

On the other hand, in Fig. 6 and Fig. 7, the situation is reverse. For solution A a = eth /la, it seems that, as the 
intensity of the pump power increases, the gain of other beams begins to saturate. Usually, when the intensity of the 
pump beams is higher than those of others, their quantum characteristics tend to vanish and the degree of quantum 
entanglement will be small. However the solution A a = e th /"f a is greater at this time, which seems to be more effective 
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FIG. 4: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7„ with 
fa = 7b = 7c = 0.03, fci = 1, fc 2 = 0.4, and e = l.le th = 
2.2eth- We use the stationary solution with A a = eth/ja- 



FIG. 5: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7 a with 
the parameters unchanged as in Fig. 4. Here we use the 
stationary solution with A a = t t h/la- 
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FIG. 6: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7„ with 
e = 2.2e th — 4.4e t h and other parameters unchanged as 
in Fig. 4. We use the stationary solution with A a = tth/'la- 
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FIG. 7: (Color online) Minima of the inequalities as a 
function of the analysis frequency normalized to 7 a with 
the parameters unchanged as in Fig. 6. Here we use the 
stationary solution with A a = t t hl la- 



va a wide range of e. 

In Fig. 8 and Fig. 9, we numerically search the minimum of V%j in all frequencies, as a function of the pump power 
parameter e, also with other parameters unchanged. In these figures, we can see clearly that the stationary solution 
A a = f~th/la is valid when e is small just like the usual cases, see Ref. 1C|. And solution A a = e th /^f a is effective 
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FIG. 8: (Color online) Maximum quadripartite entangle- 
ment as a function of the pump power parameter e with 
the parameters unchanged as in Fig. 4. Here we use the 
stationary solution with A a = eth/ja- 
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FIG. 9: (Color online) Maximum quadripartite entangle- 
ment as a function of the pump power parameter e with 
the parameters unchanged as in Fig. 4. Here we use the 
stationary solution with A a = t t h/la- 



even e is tens of times larger than e th . This directly confirms previous idea. We hope this can be considered as an 
effect of optical bistability. Although we cannot confirm the proportion of photons resulting in the stationary solution 
A a = e th /j a , it really enlarges the range of e to be effective. 



VI. CONCLUSIONS 



We propose a scheme to produce bright six-partite CV entanglement in an optical cavity operating above threshold 
in a four-level atomic system, which has been theoretically demonstrated using both VLF criterion of multipartite 
CV entanglement and the input-output relations. This result is experimentally feasible and can have a larger range 
of pump power to produce an effective multi-partite entanglement; it is much well-developed, compared with a basic 
system using the third-order nonlinearity media. Moreover, the bright six-partite entangled beams could be separated 
spatially, with narrow linewidths after generation, which have significant use in quantum network and free-space 
telcportation. 
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